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Delay probability formulas for batch input to a finite number of con- 
stant-holding-time servers are derived under the assumption of statistical 
equilibrium. The service-delay distribution (delay until a first request from 
the batch enters service) is given in terms of the roots of a transcendental 
equation, while the probability of no service-delay and the average delay are 
expressed directly in terms of the number of servers, the holding time, and 
the parameters of the input process. A numerical example with a fixed 
batch size is discussed. 

I. INTRODUCTION 

Batch arrivals constitute an important class of input processes in 
the theory of queues. The investigation of the problem of batch input 
to a group of constant-holding-time servers was motivated by the 
existence of installations with multiple Automatic Calling Units 
(ACU). Customer-based computer equipment controlling the ACU's 
is capable of originating simultaneous requests. The dial-tone markers, 
the first common control equipment in a No. 5 Crossbar central office 
to serve the requests, can be modeled as a group of constant-holding- 
time servers. 

Another example comes from an information transmission system. 
Messages containing a (small) random number of characters (a batch 
of characters) arrive according to a Poisson process and must be 
transmitted to some destination. Delayed messages are stored in a 
buffer. Since the transmission time per character is usually fixed, this 
system provides another example of the model studied. 

In Section II, the mathematical model used in this study is de- 
scribed and the input process defined ; the state equations are written 
and used to derive the generating function for the equilibrium state 
probabilities. The probability of no service-delay is found in Section 
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III, while the average delay is computed in Section IV. The service- 
delay distribution is given in terms of the roots of a transcendental 
equation in Section V. A numerical example with a fixed batch size is 
discussed in Section VI. The effect of this batching scheme on the 
average delay and the service-delay probability is examined. 



II. MATHEMATICAL MODEL 

The model studied here is that of a queuing system with a finite 
number of servers, batch arrivals, and constant holding time. The 
assumptions are 

(i) Requests arrive according to a compound Poisson process, that 
is, requests arrive in groups or batches and the instants at 
which the batches arrive constitute a Poisson process. 

(ii) There are c servers and each request has access to any one of 
them. 

(Hi) All requests have the same constant service time, t. 

(iv) The delayed batches wait until service becomes available and 
are served in order of arrival. The service discipline for requests 
within a batch is arbitrary, i.e., not specified here. 
(v) The system is in statistical equilibrium. 

Systems with simple Poisson input have been studied as early as 
1920, when A. K. Erlang 1 obtained expressions for the probability of 
delay for arbitrary values of c and the average delay for c = 1,2, and 3. 
The first complete treatment of such systems was by Pollaczek. 2,3 
Crommelin, 4 ' 5 using a method which is simpler than that of Erlang 
or Pollaczek, also derived general formulas for the probability of delay, 
the average delay, and the delay distribution. A simplified and concise 
account of Crominelin's work is given by A. Descloux, 6 who also shows 
how Pollaczek's formulas can be deduced from Crommelin's results. 
The development herein is an extension of Crommelin's results to the 
case of compound Poisson input using the simpler methods employed 
by Descloux. 

We now define the input process. Consider events which happen in 
groups rather than singly, that is, requests arriving in batches at a 
group of c servers. For k = 1, 2, • • •, let Nk{t) be a Poisson process 
with intensity A* which governs the arrival of fc-sized batches. Assume 
independence of the processes Nk(t), k = 1, 2, • • •. Let N(t) be the 
total number of requests that have arrived in the interval (0, i]. 
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Then 

00 

N(t) = £ kN k (t) (1) 

*=i 

is called a compound Poisson process (Ref. 7, page 271). 

The probability that an arriving batch is of size k is equal to X*/X, 
where 

CO 

x = £ x*. 
it— i 

From eq. (1), we see that the mean and variance of the number 
of arrivals per unit time are 

00 CO 

jui = £ k\k and n 2 = £ & 2 X*, (2) 

*=1 A- = l 

respectively. The generating function of the probability distribution 
r n (t) = P[N(t) = n], n = 0, 1, 2, • • -, is given by 

00 

v(t,z) = £ ir n (t)z n = eM'\ 
where 

00 

0(z) = £ X n2 » - X, 

71=1 

and hence the probabilities 7r„ (t) are given by 

(Xl*)*«(X 2 *) t2 ---(XnO* n 

r.(0-«-*'L- TV^ T^ ' n-0,1,2,---, (3) 

where 4„ is the class of all sets of nonnegative integers [hi, k 2 , • • • , k n \ 
such that fci -+- 2A; 2 + ■ ■ ■ + nk n = n. 

The expression given in (3) is not suited for computing x„(f). These 
probabilities are more conveniently computed from the recurrence 
relation 

t k 

T k+1 (t) = — — £ (k - j + l)X*_y +1 7r y (0, A = 0, 1, 2, • • -, 

« + 1 y=o ( 4 ) 

»o(0 = e~ X '- 

Equation (4) is easily obtained from the relation k\irk(t) = ir (fc) (i, 0), 
where the superscript denotes differentiation with respect to z. 

Special cases of the compound Poisson process are obtained by 
choosing different convergent sequences of the positive constants 
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Xi, X 2 , •••. One such sequence is obtained by setting \j = o/r'~ l , 
r > 1, j = 1, 2, •• •. This special example has become known as the 
"stuttering" Poisson process. 8 In this case, a simple expression for 
T n (t) can be obtained by noting that the generating function has a 
power series expansion in z, the coefficients of which are given in terms 
of the Laguerre polynomials L n , that is, 



e tpu) = e -\ 
since 



'H)i M -" r) Q"' 



<TZ 

1 — z/r 
It follows that 

tto(0 = e~ u 

e-x« (5) 

t»(0 = lL n (-atr) - L„_,(-^r)], n = 1, 2, • • ■. 

We will now obtain the equilibrium state equations, and find the 
probability generating function of the stationary distribution for the 
general case. Let X(t) be the number of requests (waiting or in service) 
in the system at time t. Let 

Vij {t) =P{X(t) = j\X(0) =i) 

be the transition probability functions of the process [X(t), t ^ 0}. 
It is clear that X(t) is not Markovian. If, however, we examine 

X k = X(kr), k = 0, 1, 2, ••-, (6) 

we see that this sequence is Markovian and, in fact, { X k , k = 0, 1, 2, • ■ • } 
is a homogeneous Markov chain with one-step transition probabilities 

Pij = P{X k+1 = j\X k = i\, k = 0, 1, 2, ■-., 

given by 

(ryCr), for ^ i ^ c 

Pii = hi- i+c (T), for c<ig j + c 

[0, for j + c < i. 

We will be interested in the distribution of the number of requests 
in the system encountered by an arbitrarily arriving batch (a batch 
arriving at a time point a long way from the origin, i.e., after statistical 
equilibrium has been reached). But since the instants at which the 
batches arrive constitute a Poisson process, thi6 distribution is the 
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same as the stationary distribution 

Pi = limp l7 (0, j = 0, 1,2, •••, 

i-»00 

of the process [X(t), t ^ 0} . Moreover, if this limit exists, then so does 
limPjX* = j\X Q = i] 

k -.00 

and they are equal. Consequently, the distribution of interest to us is 
given by the stationary distribution of the imbedded Markov chain 
(6). This distribution is obtained by solving the Chapman-Kolmogorov 
equations 

7>o = 7ro(r)a c 

n+c ( 7 ) 

Pn — 7T„(r)a c + £ PmTTn-m+cir), W = 1, 2, • • •, 

m = c+l 

where 

n 

«« = £ Pm- 

We will assume that h\t < c so that the stationary distribution \pj) 
exists. 

We need to solve the system (7) for the unknowns p n - To do this we 
introduce the probability generating function 

/(«) - L PnZ n - 

Multiplying both sides of (7) by z n and summing over n, we have 

/(«) = a c e*^ + £ Z p c+i r n _,(r) 2 » 

n=l j— 1 

X 00 

= a**™ + E £ p c+j ir n _j(r)z» 
;'=1 n=j 



1 

Z c 



where 



ff(2) = £ PnZ"- 

n=0 

Thus the probability generating function of the sequence p n , 
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n = 0, 1, • ■ -, is given by 

m g(z)-z<a e 
JKJ 1 -**-*« 

III. PROBABILITY OF NO SERVICE-DELAY 

We say that the service of a batch is delayed if upon its arrival all 
servers are busy. Hence, the probability of no service-delay will be 
defined by o c -i, that is, the probability of at most c — 1 servers busy. 
An explicit expression for this probability will now be found. 

We start with eq. (8). Since the p„'s are probabilities, /(z) is holo- 
morphic in \z\ ^ 1 and, therefore, the zeros of the denominator and 
numerator in \z\ ^ 1 must be the same and have the same multi- 
plicities. We will show that the denominator of (8) has c distinct 
roots in \z\ ^ 1 and that all of them, with the exception of z = 1, lie 
inside the unit circle. 

For \z\ = 1 + 5, with 5 sufficiently small and positive, we have 



I fc=i 

where mi is defined by (2). 

Since t/h < c by assumption, we have 

e r M l6+0(6 2 ) < (1 + 5) c = \Z\ C , 

and by Rouch6's theorem, the equation 

e rfi(z) _ z c - 

has exactly c roots within the region \z\ =1 + 5. Let these roots be 
denoted by z h z 2 , ■ ■ ■, z c _i, z c (= 1). Then 

(i) 1, zi, z 2 , • ■ •, z c _i are distinct. 
(it) \z„\ < 1 for n = 1, 2, ■ • • , c — 1. 

To prove (*), first note that the root z = 1 is simple because 

1 - z c e-^ (z) 
lim = Tin — c 7* ■ 

z^l z - 1 

Similarly, for any root z„ i = 1, 2, • • ■, c — 1, 

1 — Z c e~ TPill) r_l T M k "1 

lim = Zi e-^o r £ kX k z t - c • 

t _ ti z — Zi L fc=i J 
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The first two factors cannot vanish for any admissible choice of the 
root 2„ so that if z< is to be a root of second or higher order, we must 
have 

r(XiZi + 2\ 2 z 2 2 + 3X 3 2Jj + •••) = c. 

But this is not possible since 

t|XiZi + 2X 2 22 + 3X 3 z 3 + • • • I ^ rm < c, 

and the roots 1, Z\, z 2 , • • ■, z c -\ are therefore all distinct. 

To prove (w), suppose that |z„ | =1 for some n, n = 1, 2, ■ • •, c — 1, 
then |exp [t/3(z„)]| = 1 and hence the real part of T0(z n ) must be 
zero, that is, (R[/3(z n )] = 0- Hence we must have 



(R[/3(z„)] = (R 



E X fc ( 1 - z* 



= 



fc=i 



Since all terms in the sum are nonnegative, we have (R(l — z*) = 
for all k, and therefore z„ = 1, contrary to the assumption. It follows 
that |z„| < 1 for n = 1, 2, • • •, c — 1. 

Since the numerator of (8) is a polynomial of degree c, /(z) has 

the form 

A (z - \){z - zi){z - z 2 )---(z - Z c _x) 

f(z) = A — — (9) 

1 — z u e~ rfiM 

The condition /(l) = 1 determines A : 

Hit — c 
A = (10) 

(1 - 2l )(l -Z 2 )--.(l -Z C _!) 

In computing a c _i it is convenient to introduce the generating function 

F(z) = Z am n - 

Then, since a n — a„_i = p n , n = 1, 2, • • •, we have 

(1 - z)F(z) = /(z), 
or 

F(z) = ' 



1 - z 
Now, making use of (9), we obtain 

A (Z ~ Zl)(z ~ Z2)'--(Z - Z c _i) 

F(z) = —A 

1 - z'e-^ 
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The probability of no service-delay, a c -\, is given by the coefficient 
of z c_1 in the expansion of F(z) : 

"'- = ~ A = a -.0(1 -J-- -a -*<-)• (11) 

An expression for a c _i which does not involve the roots z t - is obtained 
through an application of the generalized argument-principle (Ref. 9, 
page 151). That is, suppose ^ (z) isholomorphicand0(z) is meromorphic 
on and inside the contour C. Let a*, k = 1, 2, • • •, be the zeros with 
multiplicities r k , and k , k = 1, 2, • • •, the poles with multiplicities s k 
of the function 4>(z) inside C. Then the generalized argument-principle 
states that 

— ( tW^rrdz = E nM«k) - E luKft)- 

2-iriJc <K 2 ) * * 

Taking the logarithm of eq. (11), we have 

c-l 
logflc-i = log (c - Mir) - E log (1 - Z,)- 

We will eliminate the roots Zi from the second term of the right-hand 
side of the preceding equation. Let 

(f>(z) = e T " (z) - z c , 

and note that <j>{z) has simple zeros at z = Z\, z 2 , • • •, z c -i. Choose 
}p(z) as the principal branch of log (1 - z). The generalized argument- 
principle yields 

E log (1 - Zn) = — / log (1 - z)d[log *(«)] = / 
n =i 2-jnJc 



n = l 



where C is the contour |z| = 1 — e and e(>0) is chosen so that 
z n , n = I, 2, • • •, c — 1, lie inside C but z c ( = l) is exterior to C. We 
will now show that 

co l CO 

J = log (c - 77x1) + E - E 7r;(ttr) , 
and hence 

CO 1 CO 

log a c _! = - E " E Try(nr) • (12) 

71—1 fl j=nc 

Note first that the principal branch of log (1 - z) — [z c ~ l log (1 - z)] 
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is holomorphic in | z | ^ 1 — e. Since its integral on C is zero, we have 

1 - z-<e T ^n 



1 f r i - 

J = — / log (1 - z)d log — - 
2-wi J c - 1 

Integration by parts yields 

1 r rl - 2 - c e^ (z) "| dz 

~ ~2riJc ° g l 1-2- 1 J 2 — 1 

The integrand above has a simple pole at 2 = 1, and its residue there is 
equal to log (c — hit). Choose 6(>0) such that the only zeros of 
0(2) in the disk \z\ ^ 1 + 5 are 1, z h z 2 , • • •, z c -i, and let Ci be the 
contour \z\ =1+5. Noting that the integral of log (1 — z~ l )/(z - 1) 
on Ci vanishes, we have 

J = log (c - mit) - — / log [1 - z-< e ^«>] — —■ 
2iri Joi 2 — 1 

Nowsince \z~ c e T&{ ' ) \ < 1 on &, the power series for log [1 — 2 _c e r ^<*>] 
converges uniformly on Ci, and term wise integration is allowed, so that 

1 °° 1 r zr ne 
J = log (c - mit) + — L " / « r "' (,) 7* 



= log (c - Mir) + E - 
n-i n L 



1 :£ / e Tn ^ (z) 2- nc+ ^2 ]• 



Expanding the integrand in powers of 2, and integrating term by term, 
we see that the integral is zero for all terms except one, and that there 
it is equal to 2iri times the coefficient of z~K But the coefficient of z~ x 
is exactly 7r nc _>-i(nT), so that 

Mir nc— 1 J 

J = log (C - Mir) + E ~ 1 - £ TTnc-j-linr) 

„=i n L j-a J 

co t CO 

= log (c - Mir) + L ~ L *y(nr) , 

n=l ro j'=nc 

and this is the result stated earlier. 

For the case of a simple Poisson input, we have X x = X, X, = 0, 
j = 2, 3, • • •, and (12) reduces to [Crommelin, Ref. 5, eq. (5)] 



x 1 M (Xwr)' 
log a c _i = - E ~ L — 7— 
»-in,-_ Be J! 



,-Xnr 
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Note that this expression and eq. (12) differ only in the terms it,-(»t) 
and ((XnT)Vj!)e~ XnT , which represent the same probabilities in two 
different systems: both are equal to the probability of exactly j 
arrivals in the time interval nr. As we shall see later, these probabilities 
appear again in the expressions for the average delay and the service- 
delay distribution. 

For "stuttering" Poisson input, eq. (12) reduces to 

oo e -[o-nrr/(r-l)] oo J 

loga e _i = - £ £ — [Ly(-<mrr) - Lj-ii-anrr)] 

n — 1 71 j=nc T 

where the L„(£) are Laguerre polynomials. 



IV. AVERAGE DELAY 

Under equilibrium conditions, the average delay D is equal to the 
average amount of waiting per unit of time divided by the average 
number of arrivals per unit of time. The average amount of waiting 
per unit of time is equal to 

00 

£ (n - c)p n , 

n=c+l 



so that D is given by 

1 00 1 00 

D = — £ (n — c)p n = — £ npn — T 

Ml n-c+1 Ml " = 

where mi is the average number of arrivals per unit of time denned by 
(2). An explicit expression for D can be readily obtained by noting that 

£ np n = lim [-/(z)] , |f|<l, 

n=o *-i Ldz J 

where f(z) is the generating function given by (8). Straightforward 
differentiation of (8) and application of L' Hospital's rule lead to 

D 1 c -! 1 M2T + Mit(mit - 1 ) - C(C - 1) 

— = £^ "i ~. 7 (1") 

t Mi?" i=i 1 — 2t 2mit(c — mit) 

where Z\, z 2 , • • • , z c _i are the roots denned previously. Again we wish 
to eliminate these roots. The method used in the previous section 
suggests the application of the generalized argument-principle with 
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(1 — z)~ l as yp(z), and </>(z) as before. Thus we have 

1 f 1 «'(«) 



0-1 1 

£ 

„_1 1 - Zn 2H Jcl — Z <f>(z 



d 



dz = K. 



(14) 



Noting that the integral of (1 - z)" 1 — {log |> c_1 (l — z)]} on C is 

az 

equal to (c — l)27rt (the residue at the simple pole z = is c — 1), eq. 

(14) becomes 



K - (c - 1) + — 



1 - z~ c e T ^'n 



2w; c (1 - 2) 

Integration by parts yields 

'1 — z~ c e T&< - z) ~\ dz 



k = (c- 1) — /log r — 

^ 2W7 C L 1 



- z- 1 _ (1 - z) 5 



(15) 



The integrand in (15) has a pole of second order at z = 1, and its 
residue there is found to be 



C(C — 1) — 77X1(77X1 — 1) — T/i 2 



- (c - 1) 



Consequently, 

K = - 



2(c - 77*1) 

c(c — 1) — 17x1(77x1 — 1) — 77x2 
2(c - taxi) 

1 r ri - z-<e Tf3 <*n dz 

" 2W J Cl ° g _ 1 - z- 1 J (1 -z) 2 ' 
Combining (14) and (16) and noting that the integral of 

(1 - z)- 2 log (1 - z" 1 ) 
vanishes on C\, we obtain 



(16) 



D 

T 



1 



27ri/X] 



log[l - z-«=e^(')] 



dz 



(1 - z) 2 



Recall now that \z~ e e rPlt) \ < 1 on Ci, and hence log [1 — z _e e rP(,) ] 
has a uniformly convergent power series representation in z _c e r ^ ( * ) on 
Ci, so that 



D 1 1 °° 1 T e nTt> Mzr ne 

r- / dz. 



t mit 2iri n =i n 7 c, (1 — 2) 



(17) 
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The integrand in (17) has poles at z = and z = 1 with residues 

nc—l 

£ (nc — k)ir k {riT) 

and 

n(mT — c), 

respectively, giving us the final result 

D 1 ■ 1 °° 

- - £ - £ kir nc+k (nr). (18) 

r hit „=i n fc=i 

V. SERVICE-DELAY DISTRIBUTION 

For the purpose of obtaining the service-delay distribution (delay 
until a first request from the batch enters service), let us define g m (t) 
as the probability that among the requests present at some time / , 
at most m of them are still in the system at time U + t. Considering 
the state corresponding to g mc +c-i(t), we see that, at most, mc + c — 1 
of the requests preceding the given batch will be in progress at time t 
later, and consequently, at most, c — 1 of them at time mr + t later. 
This is the condition for the service-delay d to be less than jut + t, 
or in symbols 

P\d <mr + t) = g m c+c-i(t), ^ t < t. (19) 

To determine g m (t), we introduce the generating function 

00 

G(z,t) = £ g m (t)z m . 

m— 

Upon noting that £^= ?r»_ m (0ff«(0 = a*, we have 

co oo 

G(z,t) = arww £ g m (t)z m £ ir n (t)z» 

m=0 n=0 

CO 00 

= e-«>w £ z» £ T n - m (t)g m (t) 

n=0 m=0 

oo 

= e-"»<«> £ a„z" = e- tfii '>F(z). 

71=0 

Substituting for F(z) we obtain 

1 — z e tT rpw 
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A direct expansion of the right-hand side of (20) in powers of z is 
not desirable because the coefficients of z involve sums the terms of 
which have alternate signs and, therefore, are not well suited for 
computation. To circumvent this difficulty, we first obtain a Laurent 
series expansion of G{z,t) in the annulus 1 < \z\ < |£|, where £ is 
that root of z c exp [ — t/3(z)] — 1 = which has the smallest modulus 
exceeding 1. The existence of such a root can be proved as follows. 
Since x e exp [ — t/?(x)] takes on the value 1 at x = 1 (x = real part 
of z), has a positive derivative there, and vanishes at infinity, the 
equation z c exp [ — 7-/3(2)] — 1 = has at least one root outside the 
unit circle. 

For 1 < \z\ < I £ I , the absolute value of z~ c exp [7-/3(2)] is less than 
unity. Expanding the denominator in powers of z~ c exp \_t(5{z)~] and 
the exponential function in powers of 2, and collecting like-power terms, 
we obtain, for 1 < \z\ < | £ | , 



G(z,t) = £ 



c— 1 00 

L Qn L TTcj+c+k-nKj + 1)t - <] 

11=0 y=o 



(21) 



where q n is the coefficient of z n in the polynomial 

A(z — 2i)(z — z 2 )- • • (z — 2 c _i). 

Since 2 = 1 is the only singularity of G (2, t) in \z\ < | £ | (a simple pole 
with residue —1), G(z, t) + (2 — l) -1 is holomorphic in \z\ < |£| and 
hence for \z\ < | £ | 

1 

G(z,t) H = expansion (21) + £ 2_B - 

2 - 1 n=l 

Therefore, for \z\ < 1, we must have 

00 00 

G(z,t) = L z" + expansion (21) + Er«. 

n=0 n=l 

From this equation, we obtain the service-delay distribution 

C— 1 M 

P{d < mr + t] = 1 - E ?n L 7r (m+i+2)c _ n _ 1 [(j + 1)t - f\ , 

n = ;— 

m = 0, 1, 2, •••, 0|K t. (22) 

VI. A NUMERICAL EXAMPLE 

We examine a fixed-size batching scheme which provides some 
insight into the effect of batch arrivals on the average delay and the 
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probability of service-delay. Suppose customers arrive in batches of 
size m. Then X, = for j j*m, m = m\ m and 



■7(0 = 



for j = km 



e -\mt 

fc! 

, otherwise. 



Figure 1 shows the average delay experienced for an arbitrary customer 
as a function of the occupancy p = rjui/c, for various values of m. 
We assume that the holding time is unity, and that c = 4. 
Equation (18) written in the form 



1 °° 1 

Pl n-1 n 



nc-1 



n(m - c + ce~ nXm ) + £ (nc — j)irj(n) 



was used to obtain the curves drawn in Fig. 1. We might point out 
here that the above series converges slowly when the occupancy is near 
unity. In the interest of speedy computation it may be necessary to 
solve for the roots of the denominator in eq. (8) and then use (13) 
to calculate the average delay. The same remarks apply to eq. (12) 
which is used to compute the probability of no delay. 

Because holding times are constant, several interesting phenomena 
are observed. First, if the batch size is an integer multiple of the number 
of servers, say m = kc, then the mean time until the service of an 
arriving batch (or the first customer from the batch) begins is the 
same as the average delay in a one-server system with single Poisson 
arrivals and holding time k. From the Pollaczek-Khintchine formula, 
this number is given by 

kp 



2(1 - P ) 



Hence the mean delay which an arbitrary customer experiences is the 
average of the above number and the mean delay experienced by the 
last customer in the batch to be served. Thus we have 



(m = kc). 



Note that as p -» 0, D -» (k - l)/2. 

On the other hand, if the number of servers is an integer multiple 
of the batch size, say c = jm, then the system may be viewed as a 
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p. OCCUPANCY 



Fig. 1 — Average delay in an M/D/4 queue when arrivals occur in batches of size m. 



collection of single-server systems with constant holding time and 
j-phased Erlangian input of mean offered load m/j. This can be seen 
by imagining that the sets of m servers required to serve the arriving 
batches are chosen in cyclic order. 

Figure 2 shows the probability of service-delay (the probability that 
the service of an arriving batch is delayed) as a function of the oc- 
cupancy, for various values of the batch size m and c = 4. From eq. 
(12) we obtained and used the following expression for the service- 
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Fig. 2 — Probability of service-delay in an M/D/4 queue when arrivals occur in 

batches of size m. 



delay probability : 

1-expj- L-fl- E ^(n)]l- 

Phenomena similar to those observed in Fig. 1 exist here also. For 
example, if the batch size is an integer multiple of c, say m = he, then 
the probability of service-delay is the same as the probability of delay 
in a one-server system with single Poisson arrivals, i.e., it is simply p. 

On the other hand, if c = jm, then the service-delay probability is 
the same as would be found in a system with single Poissonian arrivals 
of intensity X m and j constant-holding-time servers. 
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